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During ultrafast demagnetization of a nnagnetically ordered solid, angular monnentunn has 
to be transferred between the spins, electrons, and phonons in the system on femto- and 
picosecond timescales. Although the intrinsic spin-transfer mechanisms are intensely 
debated, additional extrinsic mechanisms arising due to nanoscale heterogeneity have only 
recently entered the discussion. Here we use femtosecond X-ray pulses from a free-electron 
laser to study thin film samples with magnetic domain patterns. We observe an infrared- 
pump-induced change of the spin structure within the domain walls on the sub-picosecond 
timescale. This domain-topography-dependent contribution connects the intrinsic 
demagnetization process in each domain with spin-transport processes across the domain 
walls, demonstrating the importance of spin-dependent electron transport between differently 
magnetized regions as an ultrafast demagnetization channel. This pathway exists independent 
from structural inhomogeneities such as chemical interfaces, and gives rise to an ultrafast 
spatially varying response to optical pump pulses. 
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Demagnetization of a ferromagnet by an optical excitation 
pulse is known to proceed on an ultrafast timescale, typi- 
cally in the femtosecond regime. Since the discovery of this 
ultrafast process by Beaurepaire et al in 1996 (ref 1), its mecha- 
nism has been intensely debated. So far, most work has focused on 
homogeneous samples to understand the intrinsic demagnetization 
mechanisms, with recent progress including the role of spin- orbit 
interaction and magnetocrystalline anisotropy^'^, the importance 
of electron-phonon-mediated spin-flip scattering^, and the inter- 
play of different demagnetization time constants at different lat- 
tice sites in ferrimagnets^. A spatial component (beyond one unit 
cell), and hence possible additional extrinsic ultrafast demag- 
netization mechanisms, has entered the discussion so far only via 
ferromagnetic-paramagnetic (FM-PM) interfaces in thin-film 
samples^'^. 

Here, we report on experiments employing magnetic thin- 
film samples exhibiting a labyrinth-like network of oppositely- 
magnetized domains with structure sizes on the nanometre 
lengthscale. As this lengthscale is also relevant for the diffusion of 
photo-excited electrons^, changes in the domain spin structure dur- 
ing femtosecond demagnetization can be expected. To achieve the 
necessary temporal and spatial resolution, our pump-probe experi- 
ments were performed at the free-electron laser (FEL) FLASH^'^ 
delivering intense, -lOOfs long soft X-ray pulses. Previous work 
combining high temporal and spatial sensitivity was either lim- 
ited in spatial resolution due to the application of a far-field optical 
probe^'^^'^^ or limited to the tens of picosecond time domain due to 
the picosecond X-ray pulses employed^^. In the present study, we 
are able to detect magnetization changes on lengthscales of a few 
nanometres with a temporal resolution of 240 fs by analysing the 
changes in the magnetic scattering signal of the X-ray probe. As a 
result, the effects of spin diffusion become experimentally acces- 
sible in a spatio-temporally matched way. We show that spin dif- 
fusion opens a topography-dependent and hence extrinsic channel 
for ultrafast demagnetization as was predicted in simulations by 
Battiato et al7. In contrast to the situation considered in refs 6,7, 
in our experiments the demagnetization effect solely relies on the 
proximity of oppositely-magnetized domains without any necessity 
for chemical or structural inhomogeneity of the sample. 

Results 

Magnetic small-angle X-ray scattering. We investigate the 
magnetic response to optical pumping for a lateral labyrinth- 
domain pattern in a thin-film Co/Pt multilayer sample as shown in 
Fig. la. Because of the strong magnetic interface anisotropy in this 
system, alternating up' and 'down magnetic domains form in a self- 
organized manner, with magnetization parallel or antiparallel to the 
sample normal. Note that the sample is laterally homogeneous with 
respect to elemental composition, that is, the boundaries between 
differently magnetized regions are purely magnetic in nature and do 
not contain spacer layers disrupting the electronic band structure. 

The sample was pumped with an infrared (IR) pulse at 800 nm 
wavelength and with 70 fs (FWHM) pulse duration. The wavelength 
of the soft X-ray probe pulse was tuned to 20.8 nm corresponding to 
the Co 3p3/2 absorption resonance, which exhibits magnetic scatter- 
ing contrast due to the X-ray magnetic circular dichroism effect^^"^^. 
Sensitivity to the spatial response during the demagnetization was 
obtained via the observation of small-angle X-ray scattering (SAXS) 
of the probe pulse. The SAXS pattern was detected on a 2D detector 
and investigated as a function of the pump-probe delay and pump- 
laser fluence, as sketched in Fig. 1 and described in the Methods. 

For SAXS from magnetic domain patterns with X-ray magnetic 
circular dichroism contrast, the scattered intensity is known to 
be proportional to the square of the local sample magnetization 
vector M projected onto the propagation vector k of the incident 
circularly polarized X-rays. Linearly polarized light, as used in our 
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Figure 1 1 Experimental set-up and data evaluation, (a) The magnetic 
multilayer sample is pumped by an optical laser pulse and probed by a 
delayed X-ray pulse produced by the FEL. The magnetic SAXS pattern 
is recorded on an IR-protected charge-coupled device (CCD) camera 
while the intense directly transmitted radiation is blocked by a beamstop. 
The inset shows a typical magnetic force micrograph of the sample in 
the probed labyrinth-domain state. The scale bar in the micrograph 
corresponds to 2jim distance, (b) By fitting the azimuthally integrated (see 
inset) SAXS intensity, the intensity maximum is determined. We observe 
a decrease in the intensity and a shift of the peak position Aqpgak, when 
comparing unpumped (blue) and pumped (green) spectra (pump fluence: 
14.2 mJcm"^, time delay: 1.3 ps). 



experiment, can be perceived as a superposition of left- and right- 
circularly polarized X-rays that have identical scattering cross- 
sections^^. The magnetic SAXS cross-section is composed of the 
structure factor reflecting the periodicity of the domain arrangement 
and of the magnetic form factor describing the mean magnetization 
profile of the domains, especially at the domain walls ^^'^^ (see 
Supplementary Discussion). As evidenced by the azimuthal symme- 
try of the SAXS ring for our magnetic system, the ensemble-average 
properties of the magnetic domain pattern are laterally isotropic, 
and we evaluate the radial scattering distribution as a function 
of the magnitude of the momentum transfer q after integration 
over the azimuthal angle 0 (Fig. lb). By fitting the intensity 
distribution, we determine the maximum intensity and hence the 
corresponding momentum transfer value for the modal value <^peak 
of the distribution. The square root of the maximum scattering 
intensity serves as a measure of the average domain magnetization 
magnitude M (ref. 14). 

Ultrafast demagnetization and spatial response. In Fig. 2a, we 
present the normalized magnetization magnitude as a function of IR 
pump fluence and pump-probe delay. These results reproduce obser- 
vations made previously for the demagnetization of single -domain 
Ni (refs 1,4), Co (ref. 4), and in particular CoPd (ref 3), thin films 
obtained by other experimental approaches: an ultrafast drop of the 
magnetization is observed with a time constant in the femtosecond 
regime, followed by a slow recovery of the magnetization on the 
picosecond timescale. The temporal resolution of 240 fs (FWHM) 
in our experiments is dominated by the arrival time jitter between 
the optical and the soft X-ray pulse. Lines represent a phenome- 
nological fit to the experimental data consisting of three different 
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Figure 2 | Ultrafast response of the magnetic system to the femtosecond 
IR pump, (a) Evolution of the intensity and (b) shift of the SAXS distribution's 
modal value in dependence on the pump-probe delay. The square root of 
the intensity serves as a measure for the domain's magnetization. Each 
data point is normalized to a SAXS measurement taken without pump. 
The typical error is indicated at a single data point for each time course. 
The time delay error is given by the time resolution of the experiment. The 
error of the normalized intensity and the peak shift were taken from the 
fit of the SAXS profile. The lines are fits to the data using the sum of three 
exponentials convolved with the time resolution of the experiment. 

exponential decays convolved with the temporal resolution, relating 
to the ultrafast demagnetization, fast recovery and long-time recov- 
ery from the quenched state^. For the initial ultrafast demagnetiza- 
tion, we achieve good fit results for a demagnetization time constant 
<300fs for all pump fluences. Subsequently, we observe a relatively 
fast but partial recovery of the magnetization within the first 2 ps 
with time constants in the range from (0.6±0.2)ps (IR pump with 
7.9mjcm"2 fluence) to (1.3±0.4)ps (14.2mjcm"2) followed by a 
slow recovery process with tens of picoseconds time constants. 

Our main experimental result is shown in Fig. 2b. Here, we 
present the relative change A^pgak/'Ipeak of the SAXS distributions 
modal value as a function of pump fluence and pump -probe delay. 
We observe a decrease of <^peak of the pumped system relative to 
the unpumped system on the femtosecond timescale, indicating 
an ultrafast spatial response of the domain system on sub-lOOnm 
lengthscales, which is observed for the first time. The tempo- 
ral behaviour is similar to that of the magnetization change, and 
we therefore fitted the spatial-response delay-time course analo- 
gously to the demagnetization. We observe an ultrafast decrease of 
^peak ^ith time constants below 300 fs towards a minimum located 
at about 0.5 ps, followed by a fast but only partial recovery of 
^peak- The recovery time constants are, within the experimental 
error, equal to the time constants observed for M. For the largest 
pump fluences around 14mjcm~^, we observe a maximum shift 
of A(^pea]^/(^pe£^]^ = (4±l)% to smaller q values, whereas we cannot 
detect any change for the lowest pump fluence of 7.9 mj cm~^. 

Discussion 

Given the inverse relationship ^ = 2K/q between the magnetic 
in-plane correlation length ^ in real space and the momentum transfer 



q in reciprocal space, such a contraction of the SAXS pattern as evi- 
denced via A^peai^/^peak by 4% would correspond to an expansion of 
the in-plane correlation length by the same ratio. In this picture, the 
average width of all up and down domains would increase by 2.8 nm 
during 0.5 ps within the entire probed area of (250x250) [im^. If 
homogeneous throughout the sample, such a hypothetical process 
must involve a rearrangement of domain boundaries over dis- 
tances corresponding to 4% of the linear field of view, that is, 10 jim 
within 0.5 ps. This would correspond to domain wall velocities on 
the order of lO^ms"^ If a domain growth process would occur in 
a hypothetic, inhomogeneous fashion such that selected domains 
were entirely annihilated or the domain width distribution was sig- 
nificantly broadened, a domain wall motion over a few nanometres 
would have to occur in the sub -picosecond regime, necessitating 
domain wall speeds on the order of 10^-10^ms~^ This speed is 
orders of magnitude higher than reported domain wall motion in 
magnetic fields or due to spin currents and exceeds the group 
velocity and even the phase velocity of magnons^^. We can thus 
rule out explanations of the observed effect based on domain wall 
motion driven by Landau-Lifshitz-Gilbert dynamics simply due to 
the fact that the resulting lateral propagation speed is too slow. On 
the same ground, a thermal expansion effect can be excluded, given 
speeds of sound below lO^ms" ^ for both Co and Pt. 

Hence, the key question is which mechanism can explain the 
observed ultrafast spatial response. We propose that spin- dependent 
transport of the optically excited electrons across the domain walls 
in our purely ferromagnetic system can provide such a mechanism. 
Experimental evidence of interlayer spin transfer has recently been 
reported for a FM-PM-FM system with ferromagnetic versus anti- 
ferromagnetic interlayer exchange coupling^. The concept of super- 
diffusive spin transport has recently been invoked in a theoretical 
study of optically pumped electrons at the FM-PM interface, with 
velocities of electrons optically excited into the sp -bands on the 
order of 10^ m s ~ ^ (ref 7). It is based on the notion that the lifetimes, 
velocities and, consequentially, the mean free paths for hot electrons 
in metals are spin dependent if the electrons move in a magnetized 
region and that secondary electrons generated in scattering events 
can contribute significantly to the resulting spin currents. Recent 
work on ballistic versus diffusive transport of spin-polarized elec- 
trons through thin epitaxial gold films is in line with this concept^^. 
As the mean free path of low-energy excited electrons in metals is on 
the order of a few nanometres, the spin-dependent electron trans- 
port does not only play a role for interlayer spin exchange between 
thin metal films but also has to be taken into account for magnetic 
systems that are laterally structured on the nanometre scale. 

To corroborate if superdiffusive spin transport can be responsi- 
ble for our experimental observations, we performed simulations 
of the lateral spin diffusion based on the work of Battiato et al7. 
We applied a Monte-Carlo simulation approach to account for the 
irregular lateral domain pattern of our sample. As a representative 
magnetization map, we used a binarized magnetic force microscopy 
image of our sample similar to the one presented in Fig. la, neglect- 
ing the width of the domain walls. As a further simplification, the 
simulation was performed on a two-dimensional lattice and spin- 
diffusion components in the direction of the sample normal were 
ignored (for simulation details, see Methods). 

Analogous to the results of the FM-PM simulation in ref 7, we 
observe a strong accumulation of minority electrons at the FM-FM 
domain boundaries. For our case of a system in which the topol- 
ogy is purely defined by magnetic order via domain walls separat- 
ing oppositely magnetized domains, this results in a softening of the 
magnetization profile at the domain walls (Fig. 3a,b). By Fourier- 
transforming the magnetization maps obtained in the simulation 
for different times after the initial optical excitation, we simulate 
the pump-probe SAXS experiment. The model reproduces the drop 
in the scattering intensity and shift in the ^peak position (Fig. 3c). 
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Figure 3 | Results from Monte-Carlo simulation of superdiffusive spin- 
transport, (a) Temporal evolution of the magnetization in a selected 
real-space region (scale bar, 200 nm) at the start of the simulation (left 
panel) and after 300fs (right panel), (b) Magnetization profiles at a 
domain boundary extracted from the magnetization maps. The blue line 
corresponds to the initial condition, the green line shows the situation after 
300fs. (c) Simulated azimuthally integrated SAXS spectrum generated 
from the initial magnetization map (blue points) and after 300 fs (green 
points). Fits (lines) were performed analogously to the ones of the 
experimental data, (d) Schematic explanation of the observed peak shift. 
Because of the softening of the magnetization profile at the domain walls 
during spin diffusion, the corresponding scattering form factor (red line) 
decreases much faster with increasing scattering vector compared with 
the initial scattering factor of a thin and steep domain wall. As a result, 
the broadened first-order magnetic structure peak in equilibrium (blue 
line) will be deformed after the pump (green line) and the position of the 
intensity maximum shifts, (e) Temporal evolution of the SAXS peak shift 
in the simulation (lines) and in the experiment (points and error bars 
taken from Fig. 2). The shaded area symbolizes the estimated confidence 
interval of the numeric results, considering the uncertainties in the input 
parameters to the simulation. Because of a convolution of the simulation 
data with the experimental time resolution, a shift is already present 
before time zero. 



As sketched in Fig. 3d, this shift is directly related to a dramatic 
change of the domains scattering form factor induced by the 
smoothing of magnetization profile at the domain walls. When 
the magnetization profile at the domain walls significantly softens, 
the corresponding scattering form factor decreases much faster with 
increasing scattering vector compared with the scattering factor of 
a thin and steep domain wall (for an explanation of the magnetic 
scattering factor, see the Supplementary Discussion). As a result, the 
structure function peak, which is strongly broadened reflecting the 
distribution of correlation lengths in the domain arrangement, will 
be deformed during spin diffusion in the domain system and the 
position of the scattering intensity maximum shifts. As this proc- 
ess is highly non-linear, a significant shift is observed only for high 
pump fluences. In our simulation model, we observe a broadening 
of the domain interface on the order of 20 nm (FWHM). 

In Fig. 3e, we compare the time evolution of the ^peak position 
shift in the simulation with our experimental data as presented in 
Fig. 2. Our model explains the sign, the timescale and — taking into 
account the large experimental error bars and the strong simplifica- 
tions of the model — the magnitude of the experimentally observed 
shift of the SAXS distributions modal value. Note that in contrast 
to the geometry in ref. 6, no interfaces between chemically different 



materials are present, and in contrast to ref. 7, it is the spin diffu- 
sion within the 2D sample plane producing the demagnetizing 
effect. In particular, no gradient in the hot- electron density due 
to the variation of the locally absorbed IR power with penetration 
depth, as always present in vertically layered systems, is required for 
this mechanism. Furthermore, the spin separation at the domain 
boundaries is, in particular, caused by the spin -dependent diffusion 
properties and cannot be explained by a simple intermixture of hot 
carriers at the domain walls, as the excitation rate in the simulation 
is independent of the electron spin. 

We have considered the ultrafast change in magnetic anisotropy 
recently reported by Boeglin et al? as an alternative mechanism 
inducing a change of the domain wall magnetization profile. The 
authors in ref. 3 find a 29% change of orbital moment, which can 
be correlated with an anisotropy change^ ^ Making conservative 
assumptions, we have carried out micromagnetic simulations with 
a 50% reduced anisotropy at t=0, taking the dynamic change in 
saturation magnetization as evident in Fig. 2 into account. Using 
a Landau-Lifshitz-Gilbert time evolution, the resulting magnetiza- 
tion profile changes by 1.2 nm after a duration of 1 ps (see Methods 
and Supplementary Fig. S2). The resulting effect is clearly too small 
to explain the observations in our experiment. 

While our simulation clearly shows the origin of the spatial 
response, the total amount of demagnetization is underestimated as 
the simulation neglects spin diffusion perpendicular to the sample 
plane^ and excludes additional spin-flip channels, for example via 
electron-phonon scattering^. More sophisticated theoretical mod- 
elling including the multilayer structure of the thin magnetic film 
and the widths of domain walls, details of the electronic band struc- 
ture and in particular spin- dependent inelastic electron scattering 
will have to be developed in the future and can be expected to lead 
to a better understanding of the effects of electronic magnetization 
transport in magnetic nanostructures. 

The observation of a spatial response of a ferromagnet— 
connected to the presence of domain walls— in response to femto- 
second optical pulses lends support to the predicted importance of 
superdiffusive spin transport in magnetically heterogeneous sys- 
tems and illustrates the generality of this concept beyond FM-PM 
interfaces''^. We note that the spatio-temporal scale of the observed 
effects is in line with the recent observation of ballistic and diffu- 
sive spin transport channels^^, when taking into account the non- 
epitaxial nature of our samples. As nanoscale heterogeneity is an 
intrinsic feature of many material systems, ultrafast effects based 
on preferential electron transport at inhomogeneities may be of 
general importance. Given the tremendous advances in the use of 
magnetic interfaces for electronic and spintronic applications in the 
past decades, one may envision the utilization of superdiffusive spin 
transport across ferromagnetic domain walls for ultrafast magneto- 
optical devices in the future. We anticipate that it will be possible to 
tailor time constants of optical demagnetization by variation of the 
domain size when approaching the spin- diffusion lengthscale. 

Methods 

Sample fabrication. The Co/Pt multilayers with a composition of Pt(50)/[Co(8)/ 
Pt(14)]i5/Pt(6) (A) were deposited on 50-nm-thick Si3N4 membranes via ion- 
assisted DC magnetron sputtering^^. By out-of-plane demagnetization before the 
experiment, the system is conditioned into a remanent labyrinth-domain state, 
which was characterized via SAXS at the BESSY II synchrotron source as well as 
by magnetic force microscopy. 

Pump-probe experiments. Experiments were carried out at beamline BL3 of the 
EEL facility ELASH in Hamburg, Germany. The X-ray probe pulse (20.84 nm) with 
a natural spectral bandwidth of about 1% was only moderately focused down to 
a spot size of 650 |im (EWHM) on the sample. We estimate an energy density of 
(2.0 ±0.6) mjcm~^ at the sample, which is below the damage threshold determined 
earlier^^. The IR (800 nm) pump pulse was delivered by a Ti:sapphire based laser 
system, which is electronically synchronized to the EEL accelerator^-^. A small 
parallax between X-ray and IR beam of 3° results in a temporal broadening of 43 fs. 
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The focus of the pump laser was set to 1.20x0. 94 mm^ (FWHM) ensuring 
nearly homogeneous pumping of the whole X-ray transparent sample area of 
250x250 |im^. The duration of the pump pulse was 70 fs (FWHM) and about 
100 fs for the probe pulse. The arrival time jitter distribution between the IR and 
the X-ray pulses had a width of approximately 200 fs (FWHM) resulting in a 
total calculated time resolution of the experiment of 240 fs. 

The magnetic SAXS experiments were performed in transmission geometry 
(Fig. la) with a charge-coupled device (CCD) area detector (2,048x2,048 pixels, 
13.5xl3.5|am^ pixel size) placed 55 mm behind the sample. The maximum 
detected scattering angle corresponds to -60 nm real-space distances, but changes 
in the ensemble-average properties of the magnetic topography can be detected 
with 3 nm accuracy as described below. The intense directly transmitted radiation 
was blocked by a central beamstop. The CCD was IR protected by a 250-nm-thick 
Al filter. Taking the filter transmission and the quantum efficiency of the detec- 
tor into account, the system had an X-ray detection efficiency of 30%. For each 
pump -probe time delay, two SAXS patterns, first without and then with pump 
pulse, were recorded from three integrated single- shot pump -probe intervals at 
1 Hz repetition rate. The pump energy density on the sample was varied between 
(7.9±l.l)mjcm-2and(14.2±2.0) mjcm'^. 

SAXS data treatment. The SAXS spectra as shown in Fig. lb were retrieved by 
azimuthal integration of the 2D patterns after background subtraction and mask- 
ing of beamstop shadow and charge scattering streaks from the membrane edges. 
All intensities were normalized to the charge- scattering signal from the membrane 
edges. Additional charge scattering was removed from the spectra by fitting a 
dependency to the intensity data at very low q where the magnetic signal is small 
compared with the charge scattering. The residual magnetic spectra were fitted 
with a split Pearson type VII distribution, which is a class of distributions that 
accounts for the skewness and kurtosis of natural distributions. In our case, the 
analytical fit function reads: 



1+ 



2 r^i,2 



(1) 



with different parameter sets (ai,mi) and (a2>^2) for the left and the right side of 
the mode ^peak- Th^ procedure results in 3 nm spatial resolution in the detection 
of changes of the in-plane correlation length. 

The delay-scan graphs in Fig. 2 were drawn from the peak intensity (Jq) values 
obtained from the measurements with IR pump (/pumped(^)) ^ ^^^^ delay 
At=t-tQ normalized to the unpumped spectra (4npumped)- Th^ delay scans were 
fitted with the sum of three different exponentials convolved with the time resolu- 
tion (a^) of the experiment (Gaussian distribution): 
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^demag> ^recovi ^recov2 denote the Strength of the demagnetization and the 
recovery process, respectively, whereas Tjemag' '^recovi ^^'^ ^recovi denote the 
time constants of these processes. H{t) is the Heaviside step function {H{t) - 1 
if t>0, H{t) = 0 elsewhere). The time zero (tg) was obtained from the fit. The 
change in the intensity peak position (^peak) is calculated analogously as 

^'2peak/^peak= ( Vmped" ^unpumped)/^unpumped and is fitted with three 
exponentials as described for the intensity change. 

Monte-Carlo simulation. The simulation starts at time zero with a number of 
0.1 excited electrons per atom corresponding to a pump fluence of approximately 
14mjcm~^ with an equal fraction of spin-up and spin-down electrons independ- 
ent of the local sample magnetization. The energy of the hot electrons is equally 
distributed over those previously unoccupied electronic states between the Fermi 
energy and the highest optically accessible state (1.55 eV above the Fermi edge). 
For each simulation run, the electron cascades initiated from 1.7x10^ primarily 
excited electrons were simulated in a lateral sample plane with a physical size of 
10xlO|im^. The magnetization maps had a spatial resolution of 4.88x4.88 nm^. 
The spin- and energy- dependent inelastic electron lifetimes were taken from 
experimental data for Cobalt^'^'^^, the Co electron velocities and all parameters for 
Pt from ah initio calculations^^'^^. The lifetimes (Tjnei) and velocities are discretized 



in channels of 0.1 eV width and averaged according to the composition of the 
multilayer. 

During the subsequent temporal evolution, we assume no interaction between 
the hot electrons and feedback of the hot electron gas on the parameters of the 
simulation is neglected as well. Without such self- consistent feedback, we expect 
our simulation to be valid only for the first hundreds of femtoseconds of the 
demagnetization process. After time zero, the electrons move in isotropically 
distributed directions on a straight-line trajectory up to an elastic or inelastic 
scattering event. The positions are updated in time steps of ty^Q = 1 fs. In case of 
inelastic scattering events, secondary electrons are generated with isotropically 
distributed, random momentum direction. The spin of the secondary electrons is 
chosen randomly independent of the scattered electron and of the local magnetiza- 
tion. The energy of the scattered electrons is calculated from classical two-particle 
collision^'^^. The elastic scattering probability (f\i) is calculated from the inelastic 
probability (tlinel = tud ^ind) as ^\ = ci^md with the constant a ranging from 

0. 5 to 2 for different MCS runs. Additional inelastic contributions, for example 
from phonon scattering are neglected. As a result of the electron-electron scatter- 
ing, the electron gas thermalizes. Every 50 fs, a magnetization map considering 
the position of the migrating electrons is saved assuming an initial magnetization 
of 1 per atom^^. 

The motion of the electrons is strongly affected by the local magnetization. 
As detailed by Battiato et al7, majority electrons possess a larger mean-free path 
length compared with minority electrons. If, however, a majority electron crosses 
a domain boundary, it will become a minority electron within the new domain 
and it is therefore trapped in the vicinity of the domain wall. As a result, these 
cascades of spin-polarized electrons within a landscape of ferromagnetic 
domains cause an enhanced demagnetization at the domain boundaries due 
to an accumulation of minority electrons. 

Micromagnetic simulations. Micromagnetic simulations of the impact of a 
change in the magnetocrystalline anisotropy were carried out using the OOMMF 
software^^, based on dynamics using the Landau-Lifshitz-Gilbert (LLC) equation. 
A 3D 300x300x36 nm^ simulation of a domain pattern was carried out with a 
cell size of lxlx4nm'^ and the resulting magnetization profile between adjacent 
domains was extracted. Values for the uniaxial anisotropy and the saturation 
magnetization have been taken from ref. 29. In the simulation, an equilibrium 
domain-wall width of 6.2 nm is yielded. Starting from the equilibrium domain 
configuration at t = 0, the anisotropy is reduced by 50% and the saturation 
magnetization by 70%. After 1 ps of LLC evolution, the domain wall broadens 
by 1.2 nm (Supplementary Fig. S2). 
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